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Abstract 

A numerical large-eddy simulation model is 
under modification and testing for application to aircraft 
wake vortices. The model, having a meteorological 
framework, permits the interaction of wake vortices with 
environments characterized by crosswind shear, stratifi- 
cation, and humidity. As part of the validation process, 
model results are compared with measured field data 
from the 1990 Idaho Falls and the 1994-1995 Memphis 
field experiments. Cases are selected that represent 
different aircraft and a cross section of meteorological 
environments. Also included is one case with wake 
vortex generation in ground effect. The model simula- 
tions are initialized with the appropriate meteorological 
conditions and a post roll-up vortex system. No am- 
bient turbulence is assumed in our initial set of experi- 
ments, although turbulence can be self generated by the 
interaction of the model wakes with the ground and 
environment. 

I. Introduction 

In response to a continuing trend for increased 
air travel that has led to more frequent delays and 
increased costs to air carriers and the traveling public, 
NASA through its Terminal Area Productivity (TAP) 
program is developing systems that will increase 
efficiency, yet ensure safety to the traveling public. A 
major element to this program is the development of an 
automated system called the Aircraft Vortex Spacing 
System (AVOSS) — which will determine safe operating 
spacings between arriving/departing aircraft as based on 
the observed/predicted weather state. 1 At the core of 
AVOSS will be a predictor algorithm that relates the 
transport and decay of aircraft wake vortices to the 
present and future weather state. The presences of 
strong crosswinds and vigorous environmental 


turbulence, for example, may lessen the threat of a 
dangerous wake vortex encounter, and thus allow a 
tightening up of wake vortex separation standards now 
in effect during instrument meteorological conditions. 

NASA-Langley theoretical modelling efforts are 
expected to play a major role in the development of 
these algorithms. A numerical large-eddy simulation 
model called the Terminal Area Simulation System 2 
(TASS) has been modified for application to aircraft 
wake vortex simulations. The TASS model, which is 
described in section 2, is believed to have an advantage 
over other wake vortex models 3 ’ 4-5 ’ 6 ’ 7-8 in that it has a 
meteorological reference frame, compressible non- 
Boussinesq equation set, subgrid turbulence closure, a 
formulation for ground-friction, realistic boundary 
conditions, an option for either two or three spatial 
dimensions, and accurate — yet computationally efficient 
— numerical approximations. In spite of the model's 
sophistication, it has relatively fast execution times, 
which is an important feature when taking on the 
ambitious task of simulating three-dimensional, time- 
dependent aircraft wakes. The TASS model is capable 
of simulating post roll-up wake vortices in both two and 
three dimensions, for a wide range of atmospheric 
conditions that include: vertical wind shear, 

stratification, atmospheric boundary layer turbulence, 
fog, and precipitation. Initial emphasis is being placed 
on validation. Once this has been accomplished with 
sufficient confidence, then TASS model output will be 
used to develop parametric relationships for vortex 
transport and decay as related to different aircraft types 
(e.g., B-767, DC-10, etc.) and meteorological conditions 
(windshear, stratification, ambient turbulence, humidity 
and precipitation). The model can also provide high 
spatial resolution fields of wind, pressure, temperature, 
and humidity, for analyzing and understand the structure 
of wake vortices. The generated high-resolution fields 
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will be valuable in the characterization and quanti- 
fication of wakes generated in various meteorological 
scenarios, as well as for providing realistic data sets for 
wake vortex sensor development. 

The purpose of this paper is to give an over- 
view of the TASS model and to present results of its 
application to wake vortex research. Two-dimensional 
simulations of cases from the 1990 FA A Idaho Falls 
and from the Memphis 1994 and 1995 NASA/MIT- 
Lincoln-Laboratory wake-vortex field deployments are 
presented. Results from the computational simulations, 
including comparisons with measurements, are presented 
for the Idaho Falls cases in section 4 and for the 
Memphis cases in section 5. Details of the TASS 
model are described in the next section and the initial- 
ization procedure is reported in section 3. 

II. Model Description 

The TASS model is a multi-dimensional, large- 
eddy code that has been used to study a variety of 
meteorological phenomenon including convective local 
storms, 9 microburst/windshear, 10111213 hailstorms, 14 
tornadic thunderstorms, 15 nuclear cloud rise, 16 and 
atmospheric boundary layer turbulence. 17 ' 18 The TASS 
model has been validated with observed data for a 
range of micro- and cloud-scale phenomena, and was 
highly successful in the completed NASA-FAA wind- 
shear program. Recent modifications to the initial 
conditions allow for its application to aircraft wake 
vortex phenomenon. Thus the model is not merely a 
"wake vortex model," but can be generically applied to 
a diversity of phenomena. The model includes parame- 
terizations for ground stresses that are a function of 
surface roughness, allowing for "in-ground effect" wake 
simulations with realistic ground interactions. 

The TASS model consists of 12 prognostic 
equations: three equations for momentum, one equation 
each for pressure deviation and potential temperature, 
six coupled equations for continuity of water substance 
(water vapor, cloud droplet water, cloud ice crystals, 
rain, snow and hail) and a prognostic equation for a 
massless tracer. The non-Boussinesq formulation used 
in TASS was initially developed to study tornadoes with 
application to similar intense vortices. 19 This 
formulation was subsequently extended to a com- 
pressible time-split formulation, and parameterizations 
for numerous microphysical interactions were included 
as an option. 2 Salient characteristics of the TASS 


Table 1. Salient Characteristics of TASS 6.x 


Primitive equation / non-Boussinesq equation set 
Time-dependent, nonhydrostatic, compressible 

Meteorological framework with option for either 
three-dimensional or two-dimensional simulations 

Liquid and ice phase microphysics -- can simulate 
growth processes for cloud droplets, cloud ice 
crystals, rain, snow and hail. Accomplished with 
large set of microphysical-parameterization models 

Inverse-exponential size distributions for precipitat- 
ing hydrometeors 

Large Eddy Simulation model with Ist-order sub- 
grid scale turbulence closure - scales of turbulence 
larger than grid volume are resolved in the 
simulation 

Ground stress based on Monin-Obukhov Similarity 
theory 

Choice of lateral boundaries: may be either open, 
mirror, or periodic - open condition utilizes mass- 
conservative, nonreflective radiation boundary 
scheme 

Option for nonstationary domain - movable, 
storm/vortex centering mesh 

Explicitnumerical schemes, quadratic conservative, 
time-split compressible- accurate and highly 
efficient, almost no numerical diffusion 

Arakawa C-grid staggered mesh, and vertical 
coordinate stretching allowed 

Filter and sponge applied to top four rows so as to 
diminish gravity wave reflection at top boundary 

Ambient atmospheric conditions initialized with 
vertical profile of pressure or altitude, temperature, 
dew point, and wind velocity 

Model applicable to meso-y and micro scale 
atmospheric phenomenon. Initialization modules for 
simulation of convective storms, microbursts, 
atmospheric boundary layers, and aircraft wake 
vortices 


model are listed in table 1. 
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The TASS equation set in standard tensor 
notation is as follows: 


Conservation of Scalar Variables (e.g., water vapor, 
cloud droplet water, etc.) 


Momentum: 


du, 

dt 


H dp du i u j 

— = + u . 

Po (h i dx j 

-2Q.J (u k ~u 0k )E ijk 


du 

■+ g(H -1)5 j3 
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— = + — + + Sq 
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where Sq represents contribution from source terms. 
For precipitating variables (such as rain and snow), an 
additional vertical flux term is added to the above 
equation to account for fall out. 


Buoyancy Term: 
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Pressure Deviation: 
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with the Potential Temperature being defined as: 


0 


Poo \ 


T (—) 
P 


In the above equations, u i is the tensor component of 
velocity, t is time, p is deviation from atmospheric 
pressure P, T is atmospheric temperature, p is the air 
density, Q. is the earth's angular velocity, C p and C v are 
the specific heats of air at constant pressure and volume, 
g is the earth's gravitational acceleration, R d is the gas 
constant for dry air, P ()() is a constant equivalent to 1000 
millibars (10 5 pascals) of pressure, Q v is the mixing 
ratio for water vapor, and Q T is the sum of the mixing 
ratios for liquid and ice water substance, L v is the latent 
heat for vaporization for water, L f is the latent heat for 
fusion for water, L s is the latent heat for sublimation 
for water, and s v , jy , and s s are the respective water 
substance source terms. Environmental state variables, 
e.g., u ok , Q vo , p o , and Q a , are defined from the 
initial input sounding and are functions of height only. 


The dependent variables in TASS are treated as 
averages over the grid volumes — giving rise to subgrid 
Reynold's stress and subgrid eddy transport terms. 
From first-order closure theory, the subgrid scalar 
covariances and subgrid stress are respectively: 

5/Q) =p d K h -g 


x ij 


du . 3 u 

■ Po 

OX; OX; 


2 du k 5 

3 3.r ij 


J = P 0 K m If, 


A modified Smagorinsky first-order closure is used for 
the subgrid eddy viscosity: 


K M =’s 


N 


du du du 
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du 

dx 


-) 2 



The subgrid eddy viscosity is modified by stratification 
through the Flux Richardson Number, If- , which is 
approximated from the Richardson number, Ri, as: 


R f = Min[ 0.9998, Ri K H /K M ], 

where the subgrid eddy viscosity for heat is specified as 
K n = 3 K m . The Richardson number for moist, but 
unsaturated air is 20 

1 dP 30(1 + 0.61g v - Q t ) 

P» 9 dx k dx k 


For saturated air (i.e. 100% relative humidity) the 
Richardson number is modified from above to account 
for the moist adiabatic ascent and descent of saturated 
air parcels. Therefore, (l/Q)dQ/dz is replaced with: 
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J_5e _ g 

e dz C p T 


where R v is the gas constant for water vapor, e = 
R/R v , z is the vertical coordinate, and e v is the vapor 
pressure for liquid water. The subgrid length scale, l s , 
is determined from the grid volume and is matched to 
the appropriate length scale close to the ground where 
the flow is under-resolved. That is: 


£ Vv 

R d T 


C p R r T 2 


aA 

aA[l + ( aA/kz ) m_1 ] 

1 + (aA IkzT 
kz 


z > aA Ik 
a Aik > z > Az/2 
z < Az/2 


where k is von Karman's constant, and where m and a 
are invariant constants with values defined as m = 2.5 
and a = 0.155. The filter width is 


A = [2Ar 2Av 2Az ] 1 /3 

where Ax, Ay, and Az are the numerical grid sizes in the 
respective x, y, z direction. 


generating aircraft. In the event of strong crosswinds, 
wide domains are not needed since the numerical grid 
translates with the horizontal advection of the wake 
vortices. The top of the domain is chosen at least 50 m 
above the altitude of the generating aircraft. In the 
following cases, the horizontal domain size is between 
150 m and 300 m wide and the vertical depth is be- 
tween 100 m and 250 m. The simulations use the open 
radiation boundary condition option at the two lateral 
boundaries, allowing minimal interference and distortion 
to the interior flow. 

Ambient Conditions 

Each case is initialized with its respective 
vertical distribution of observed temperature, dewpoint, 
and wind velocity, that represents the air mass surround- 
ing the wake. Horizontal variations in the ambient 
conditions are not considered. 

Vortex Initialization 

As for the initial wake field, the two-dimen- 
sional simulations are initialized with a simple vortex 
system that is representative of the post roll-up, wake 
velocity field. The vortex system is initialized with the 
superposition of two counter-rotating vortices, with the 
velocity field for each vortex according to the Burnham- 
Hallock model 21 as: 


The ground boundary is impermeable (except 
to precipitation) with nonslip velocity specifications. 
Surface stresses are computed locally from the near- 
ground wind velocity, Richardson number, and surface 
roughness height. 

An abbreviated discussion of the numerical 
approximations for the TASS equation set can be found 
in the Appendix. 


where V is the vortex tangential velocity, r the radius, 
r c is the core radius (i.e. radius of peak tangential 
velocity), and Y v the circulation at r > > r . Appropri- 
ate vortex image conditions are applied to the initial 
wake field to ensure consistency and mass continuity at 
the model boundaries. 


III. Initial Conditions 

Grid/Domain Configuration 

For the simulations presented in this paper we 
have chosen the 2-D option for all simulations with Ac 
= Ay = Az , and have rotated the coordinate system 
such that y is along the aircraft flight path and x in the 
cross track direction; z represents distance above the 
ground. The simulations are time integrated within the 
two-dimensional x-z plane. In the cases presented here, 
a constant grid size of 3/4 to 1 meter is used. 

The computational domain size for each case 
simulation is governed by the size and height of the 


The parameters that govern the vortex initial- 
ization are: vortex separation, core size, height, and 
circulation. The height of the vortex system, Z ( , is 
determined from the observed height of the generating 
aircraft. The vortex core radius is assumed to be 5% of 
the generating aircraft's span. The remaining parameters 
are aircraft dependent and can be obtained from aerody- 
namic theory with the assumption of an elliptically- 
loaded wing 22 as: 

S = —B = 4Mg 

4 7ipC fl S 
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where S is the vortex separation, B is the span of the 
generating aircraft, V u is the true air speed, and M is the 
mass of the generating aircraft. 

This initialization procedure does not account 
for the aircraft configuration (i.e. flap settings, landing 
gear, inclined flight path); however, results from this 
procedure were found to give good comparisons with 
observations regardless of flight configuration. 

Turbulence Initialization 

All of the following cases assume that there is 
no preexisting ambient turbulence. However, turbulence 
(and the effects of turbulence) can develop in the 
simulations from the interaction of the wakes with either 
the ground or ambient environment. This assumption of 
no initial turbulence is acceptable for stably-stratified 
environments, since the affect of vortex generated 
turbulence likely overwhelms the affect due to any 
preexisting turbulence. Currently, an initialization for 
ambient turbulence is being developed, and future 
studies will evaluate its effect in both two- and three- 
dimensional wake simulations. 

Execution Time 

The TASS code can be run on any modern 
computer or workstation that is equipped with a fortran 
compiler. However, the domain size and grid resolution 
will be limited by the host computer's memory and 
speed of execution. A typical two-dimensional wake 
vortex simulation with a 1-m grid resolution takes about 
20 min of CPU time on a Cray C-90 supercomputer. 
For a fixed domain, the execution time is inversely 
proportional to the cube of the grid size. So a doubling 
(halving) of the resolution can increase (decrease) the 
execution time by a factor of eight. 

IV. Results — Idaho Falls 

The first set of TASS validation experiments 
are from the 1990 Idaho Falls, Idaho (IDF) field 
experiment. The principal aim of this FAA sponsored 
field study was to evaluate Boeing 757 and 767 aircraft 
wakes in a week-long series of tower fly-bys. An 
extensive data set was generated by a DOT/NO AA team 
containing wake vortex measurements and atmospheric 
data. 23 ’ 24 Vortex data was collected from a series of 
low-level flights in various flight configurations by a 
Faser Doppler Velocimeter (FDV), an array of Mono- 
static Acoustic Vortex Sensing System (MAVSS) 
sensors, and a 200 ft (61 m) instrumented tower. 


Atmospheric data were measured by the tower and a 
tethersonde. The tests were conducted in September 
during mornings and extending into early afternoons. 
The prevalent atmospheric conditions were dry and 
stable with local drainage flows — transitioning as the 
morning progressed — to an unstable boundary layer that 
grew in depth due to solar heating. 

From this data set, seven cases were chosen 
that represent a range of atmospheric conditions and 
aircraft types (table 2). Characterization of the meteorol- 
ogical conditions are represented by several parameters 
(also listed in table 2), which are: the vertical change in 
potential temperature A0/Az (negative for unstable 
stratification, positive for stable), mean crosswind shear 
AU/Az , crosswind at flight level U(Zj) , and a Bulk 
Richardson number based on the crosswind flow, i.e., 

Ri b - i - A9/A2 
9 (AU/Az) 2 

where the mean gradients in the listed parameters are 
taken between 5 m AGF and 15 m above the flight 
path. 

Model validation experiments were conducted 
for each of the seven cases, with the input for ambient 
conditions constructed from the tower and tethersonde 
data. The tower provided excellent ambient wind and 
temperature information up to 61 m. At higher eleva- 
tion the tethersonde data was used, which unfortunately, 
was of lesser quality, being subject to errors from very- 
short averaging times and bobbing of the balloon. 
Values for the initial vortex parameters (also included in 
table 2) were determined from aircraft dependent 
specifications and the reported weight, height, and speed 
of the generating aircraft. Comparison of the simula- 
tions with the measured data was generally very good 
with computed vortex trajectories having the best 
agreement. Good agreement was obtained independent 
of flight configuration and aircraft type. Comparison of 
model circulation values tended to upper-bound the 
observations, with the measurements indicating greater 
decay rates after one or two minutes. This was better 
than expected since the two-dimensional simulations do 
not allow vortex stretching and do not permit develop- 
ment of 3-dimensional decay processes such as vortex 
bursting and Crow instability. Three of the case 
simulations (listed in bold in table 2) are presented 
below. 
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B-757 Run-9 

This particular case was selected by the FAA 
as a standard against which all wake-vortex modelers 
could compare. 25 The field observations indicated a 
strong and persistent upwind vortex from a B-757 fly-by 
at 0818 MST on 25 September. The B-757 was in a 
landing configuration and on a level flight path at 70 m 
AGL. The ambient atmosphere for the early morning 
flight was characterized by significant crosswinds (5.8 
m/s at flight level). The lapse rate for ambient tempera- 
ture was stable, except below 20 m AGL due to the 
influence of the ground beginning to heat up. A strong 
crosswind shear was especially noted between 20 and 
40 m AGL. 

As shown in Fig. 1, comparisons between field 
data and TASS model results are extremely good. Both 
observations and model results show that the vortices 
initially descend due to the mutual interaction of their 
velocity fields. Flowever, the vertical descent of the 
downstream (port) vortex is suppressed as it encounters 
moderate crosswind shear of opposite vorticity; and it 
eventually ascends upward with increasing lateral 
separation from the upstream (starboard) vortex. The 
upstream vortex descends to and remains near the 
ground, where it translates at relatively slower speeds 
due to the presence of weaker crosswinds. The model 
results show several minor bounces of the upstream 
vortex, resulting from the vortex interaction with the 
ground and production of secondary vortices. 

The wind vector and potential temperature field 
at 90 sec are shown in Figs. 2 and 3, with Fig. 2 also 
showing the observed positions of the wake vortices. 
Note that potentially warmer air is carried downward 
with the downwash between the vortices, and potentially 
cooler air is transported around the periphery of the 
vortices. [Potential temperature is conserved during dry 
processes, and thus acts as a fluid tracer.] Also evident 
in Fig. 3, is the presence of complex interactions 
between the ground and the upstream vortex, including 
the presence of an obvious secondary vortex. 

Average Circulation 

In order to characterize the strength of wake 
vortices, previous field investigators have defined a 
parameter called average circulation. 21 According to 
reference 21, this parameter is preferable since it relates 
to the rolling moment of an encountering aircraft and it 
provides a more stable measurement than circulation. 
The average circulation is defined as 



where the circulation is defined as 
F =£ V-dl 

or equivalently, 

r ' -jit m 

where a and b are the radii of the averaging interval, 
and £ is the axial component of vorticity. Values for 
these parameters are easily determined from simulation 
data, but less so easily from observations. 

Fig. 1 shows a comparison of the 10 m average 
circulation (i.e. the average circulation from r — 0 to 
r = 10 m) for TASS data, as well as that deduced from 
the tower, LDV and MAVSS measurements. The 
modeled 10 m average circulations of the two vortices 
are about equal in magnitude, and bound the estimated 
circulation from field measurements. The average 
circulation estimated from the MAVSS sensor seems too 
low, especially when compared with the LDV and tower 
data. The MAVSS estimates were probably low due to 
its poor resolution of the vortex core. 26 

Fig. 4 shows average circulation vs averaging 
radius at two different times from Run 9. The com- 
parison between TASS and LDV data are quite good, 
even near the core region of the vortex. 

Comparisons in Fig. 5 from TASS data show 
that the magnitude of the 5-15 m average circulation is 
closer to that of the initial circulation 
[T(r = 0, r >> r ( )\. Also apparent from the normalized 
curves in Fig. 5, the 10 m average circulation decays 
faster than the 5-15 m average circulation. 

The purpose of the above discussion on circula- 
tion is not to advocate one parameter over the other, but 
to show that they can differ in magnitude and in rates of 
decay. When validating model simulations and defining 
separation standards based on vortex decay rates, one 
should not assume, for example, that the 10 m average 
circulation is (and decays) approximately the same as 
circulation. 
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IDF B-767. Runs 23 and 3 1 

Contrasting with Run 9 are two cases from B- 
767 fly-bys that were flown in environments of weaker 
crosswind shear. In Run 23 the aircraft was in "take- 
off" configuration within a stably-stratified environment. 
While in Run 3 1 the aircraft was in "landing configura- 
tion" within an environment of unstable stratification. 
Comparison between the TASS simulations and field 
data for these cases are shown in Figs. 6 and 7. 

For run 23 (Fig. 6) there are no observations of 
the downstream vortex after 60 sec. In the model 
simulation after 80 sec, the downstream vortex dimin- 
ished in size and intensity, as it was coalesced within 
the larger-scale circulation of the downstream vortex. 

In run 3 1 (Fig. 7), the downstream vortex was 
no longer tracked after 40 sec, although the simulated 
downstream vortex remained intense. 

In both of these events it is not known whether 
the downstream vortex suddenly decayed or was lost by 
the LDV sensor. In each of the events, both upstream 
and downstream vortices did encounter the tower, at 
positions and time intervals in agreement with the model 
simulation. 

V. Results -- Memphis 

The last set of TASS experiments are from the 
Memphis field experiments. The NASA-Langley 
sponsored field study was conducted by MIT Lincoln 
Laboratories, in November and December of 1994 and 
again in August 1995. 27 28 The purpose of the study 
was to gather meteorological, wake vortex, and aircraft 
data at an operational airport, for use in validation of 
wake vortex models and for direct use in the develop- 
ment of the AVOSS prediction system. Aircraft wake 
vortex measurements were obtained with a 10.6 micron 
continuous wave laser with real time identification and 
tracking algorithms. 29 Meteorological data included: a 
150 ft (46 m) tower with sensors for wind velocity, 
temperature, and humidity; a radar profiler and acoustic 
sodar for measuring winds aloft; a radio acoustic 
sounding system (R ASS ) for providing temperature data; 
and frequent rawinsonde balloon launches, for obtaining 
the vertical profiles of wind velocity and temperature. 
Aircraft data, such as weight, type, and airspeed, as 
well as aircraft beacon data were collected. Lidar 
measurements were limited to arriving aircraft, with 


some of the 1995 cases set-up for measuring wakes 
generated in ground effect. 

Simulation results and lidar data are shown in 
Figs. 8 and 9 for four of the cases. The selected cases, 
which are summarized in table 3, represent a cross 
section of aircraft commonly encountered at 
international airports. One of the cases (#1475) is for a 
MD-11 wake generated in ground effect. 

The meteorological conditions for the selected 
cases are characterized as having light crosswinds with 
weak vertical shear. Since the shear was light, the 
vortex pairs descend to the ground with little increase in 
separation. The crosswinds for cases #1254 and #1475 
are especially light, resulting in little lateral movement 
of the wakes (Fig. 9). 

Comparisons between the modeled wake 
trajectories and lidar measurements are excellent, even 
for case #1475, which is the MD-11 wake generated in 
ground effect. Both the TASS simulation and lidar 
measurement show the MD-1 1 wake slowly rising from 
its generation altitude of 17.5 m to an elevation of about 
40 m (Fig. 8), with the starboard vortex remaining in 
the vicinity of the flight path (Fig. 9). Unfortunately, 
lidar measurements of the port vortex were unavailable 
for this case and two of the other Memphis cases. 

VI. Summary 

Two-dimensional TASS wake vortex simula- 
tions were compared with field measurements for a 
range of meteorological conditions and aircraft types. 
Successful validation of two-dimensional TASS results 
were achieved with the initialization of the appropriate 
observed meteorological conditions and a post roll-up 
vortex system. Excellent agreement is obtained between 
TASS predicted trajectories and measurements. Good 
agreement was obtained whether the generating aircraft 
was in landing or take-off configuration. Results from 
the simulations indicate great sensitivity to meteorology, 
especially vertical wind shear. A vortex encountering a 
shear with opposite sense vorticity can be deflected 
upwards. Circulation values predicted by TASS are 
upper-bound for the observed values. The field data 
indicated a greater decay rate once the wake vortex was 
one to two minutes old. However, the future inclusion 
of initial turbulence fields and the extension to three 
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dimensions should lead to a more realistic treatment of 
vortex decay. 


Following successful validation, TASS para- 
metric runs are to be used to quantify the effect of 
weather conditions on wake-vortex transport and decay, 
and provide data useful for development of predictor 
algorithms for aircraft spacing. 


Table 2. Idaho Falls Validation Cases. Only those cases in bold are presented in this paper. 
[A0/Az is the change in potential temperature with height, AU/Az is the change in crosswind 
with height, and U(Zj) is the crosswind speed at the generation height] 


IDF 
Run # 
& date 

Aircraft 
& Con- 
figuration 

Meteor- 

ological 

Conditions 

Initial Vortex 
Parameters 

Environmental Parameters | 

(nr/s) 

5 

(m) 

Z/ 

(m) 

Bulk 

Rich- 

ardson 

A6/A- 
(°C per 
100 m) 

AU/Az 
( 10" 2 s" 1 ) 

U(Zj) 

(m/s) 

#6 

9/23 

727-222 

landing 

stable / 
moderate 
shear 

320 

26 

79 

1 

8.5 

5.2 

5.0 

# 22 
9/23 

727-222 

landing 

near- 
neutral / 
moderate 
shear 

300 

26 

76 

.04 

.1 

3.0 

5.4 

# 9 
9/25 

757-200 

landing 

stable / 
moderate 
shear 

365 

30 

70 

.5 

3 

■ 

5.8 

# 30 
9/25 

757-200 

landing 

unstable / 
low shear 

360 

30 

70 

■ 

-0.9 

0.7 

2.8 

# 45 
9/26 

757-200 

takeoff 

stable / 
low shear 

355 

30 

79 

1.9 

5.9 

3.2 

2.8 

# 7 
9/29 

767-200 

landing 

stable / 
moderate 
shear 

375 

38 

70 

0.6 
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Figure 1 . Comparison of TASS results with field data for the Idaho Falls B-757 Run 9 case. The figures (left to 
right) are: altitude of the vortex track with time, lateral position (relative to tower) vs time, altitude vs lateral 
position, and 10 m circulation vs time. Upstream positions indicated by solid lines (simulation) and filled symbols 
(measurements), downstream positions indicated by dashed lines (simulation) and open symbols (measurements). 
Right triangles, squares, and circles represent LDV, MAVSS, and tower data, respectively. The average circulation 
estimates for the upstream and downstream vortex from MAVSS are: M and m , respectively. [Field data courtesy 
of Volpe National Transportation System.] 
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Figure 2. Cross section of simulated wind vector field for IDF B-757 Run 9, at 90 sec after vortex generation. 
Aircraft track at: x = 0 and z = 70 m , tower position at: x = -95 m, both being to the right of the domain. Wind 
vectors are relative to movement of vortex pair. Observed vortex locations are denoted by filled circles. 



Figure 3. Same as Fig. 2, but for simulated potential temperature. 
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Figure 4. Average circulation vs averaging radius for IDF B-757 Run 9, at 30 sec and 60 sec. Simulated 


data represented by solid line and LDV data by solid circles. 




Figure 5. Simulated average circulation vs time for IDF B-757 Run 9. Comparison of 10 m average 
circulation and 5-15 m average circulation. Right figure depicts normalized circulation history. 
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Figure 6. Same as Fig. 1, but for IDF B-767 Run 23. 
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Altitude vs Time Lateral Displacement vs Time 






Figure 7. Same as Fig. 1. but for IDF B-767 Run 31. 


Table 3. Memphis Validation Cases (symbols as defined in Table 2) 
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2/ 
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Rich- 

ardson 

A6/A z 
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# 560 
DC-9 
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unstable/ 
weak shear 
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-800 
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~0 

1.9 
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DC- 10 
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480 

37 
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# 1254 
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8/16/95 
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stable/ weak 
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26 
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0.5 

# 1475 
MD-11 

8/25/95 

0344 

stable/ 
weak shear 

560 
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1.6 
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DC-9, 1 4 Dec 1 994, 0113 UTC (#560) DC-10, 1 4 Dec 1 994, 0630 UTC (#586) 




B727-1 00, 16 Aug 1995, 0319 UTC (#1254) MD-11, 25 Aug 1995, 0344 UTC (#1475) 




Figure 8. Comparison of TASS results with field data for the four Memphis cases. The figures are altitude of the 
vortex track vs time with simulated starboard and port portions indicated by solid and dashed lines, respectively. 
Right triangles and squares represent starboard and port lidar measurements, respectively. 
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Figure 9. Same as Fig. 8, but for lateral position of the vortex track vs time. 


Table A-l. TASS 6.x Numerics 


Prognostic Variable 

Time Derivative 

Space Derivatives 

Momentum and 
Pressure 

Time-split 30 — small time step 
for acoustic terms, 

2nd-order Adams-Bashforth: 
both large and small time steps 

Centered, Quadratic-Conservative Dif- 
ferences — with 4th Order Accuracy for 
Convective terms, remaining terms 
2nd Order Accuracy 

Potential Temperature, 
Water Substance, etc. 

Third-Order time/space with Upstream-Biased Quadratic Interpolation 31 


Appendix 

TASS Numerics 

Factors that influence the choice of numerical 
schemes should be: accuracy, economy, and long-term 
stability; truncation error should not allow the artificial 
accumulation (or dissipation) of mass, momentum. 


vorticity, and potential and kinetic energy. Numerical 
dissipation and/or moderate applications of explicit 
filtering should be avoided in computational simulations 
of wake vortex phenomenon, since it may result in the 
artificial dilation of the vortex core and in exaggerated 
vortex decay rates. The time-derivative approximations 
in TASS (see Table A-l ) are explicit (for computational 
efficiency) with minimal numerical dissipation. The 
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spatial schemes in TASS are — Quadratic Conservative 
Schemes 32 — which obey certain integral constraints 
on quadratic quantities, such as kinetic energy. The 
convective derivatives in the equation of motion are 
approximated on an Arakawa C-grid 33 (staggered in 
space) with 4th-order centered finite differences as 
formulated by Wilhelmson and Chen 34 . 

Using standard finite difference operator 
notation as 


8 „vG(-u) - 


i Ax 


[<2Cv. 


i Ax . 


ru , 

Q(xr — — )] 


and the finite difference form for the flux component 
for the w velocity is: 

(uw) x + (vw) y +(ww) z “ 1^-5 x (u z w x ) 

- Is r (u z w 3x ) - Is, (u 3z w x ) + i°b V (v z w y ) 

8 A 8 3a 8 y 

- is (v z w 3y ) - is, (T 3 -' w y ) + i^S (w z w z ) 

8 A 8 3 > A 8 z 

— -!-5-(w 3 '' w 3 ') - -i-Sg^w'' w z ) 

8 8 

The divergence term in the w-equation of motion in 
finite difference form is: 


m^) nAx - -UQ(X, + + Q(x, - ^1)] 

2 2 2 


nAx , 


nAx . 


u[u x +v y + wj - i([_(5 v i( + S y v +5 r w) 

1 


_(5 3v i< + S 3y v +5 3 ; w)] 


the finite difference form for the flux component for the 
u velocity is: 

( uu ) +(vm) + (w«). - -4 -S..(m a u x ) 

a y 4 Q a 

- ^,(« 3A « 3A ) - «*) + ^,(v A u y ) 

- hp x u 3y ) - h,p 3x u y ) + }^ z (w x u z ) 

- _W u 3z ) - -^5 3z (w 3 ' c Tr) 


and similarly for v an w. 


Similarly, the finite difference form for the flux compo- 
nent for the v velocity is: 


(uv) x + (vv) + (wv) z - ^ x (u y v x ) 

o 

- v 3A ) - ^3 > 3v v a ) + ^ y ( V y v y ) 

- ^ y (v 3y V 3y ) - v y ) + ] h z (w y v z ) 

O' O' o 
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